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MAIN BODY 

Like a quantum computer designed for a particular 
class of problems, a quantum simulator enables quan- 
titative modeling of quantum systems that is computa- 
tionally intractable with a classical computer. Quan- 
tum simulations of quantum many-body systems have 
been performed using ultracold atoms^ and trapped ions^ 
among other systems. Superconducting circuits have re- 
cently been investigated as an alternative system in which 
microwave photons confined to a lattice of coupled res- 
onators act as the particles under study with qubits cou- 
pled to the resonators producing effective photon-photon 
interactions". Such a system promises insight into the 
noncquilibrium physics of interacting bosons but new 
tools are needed to understand this complex behavior. 
Here we demonstrate the operation of a scanning trans- 
mon qubit and propose its use as a local probe of photon 
number within a superconducting resonator lattice. We 
map the coupling strength of the qubit to a resonator 
on a separate chip and show that the system reaches the 
strong coupling regime" " over a wide scanning area. 

Over the past decade, the study of quantum physics us- 
ing superconducting circuits has seen rapid advances in 
sample design and measurement techniques"". A great 
strength of superconducting qubits compared to other 
promising candidates is that they are fabricated using 
standard lithography procedures which allow fine tun- 
ing of qubit properties and make scaling up the fabrica- 
tion to devices with many qubits straightforward. Circuit 
quantum electrodynamics (CQED) is an active branch of 
quantum physics research in which one or more qubits are 
strongly coupled to a superconducting coplanar waveg- 
uide resonator (CPWR) which is used to control and 
readout the state of the qubits"". A prerequisite for 
most interesting CQED applications is that the system 
reach the strong coupling regime in which the rate at 
which the qubit and the resonator exchange an excita- 
tion exceeds the excitation decay rate. 

In addition to the CQED architecture's promise as 
a quantum computing platform, recent theoretical work 
has focused on using a CQED lattice, a network of cou- 
pled resonators each coupled to its own qubit, as a non- 
equilibrium quantum simulator". One particularly in- 
teresting prediction for CQED lattice systems is a cross- 
over from a superfluid-like state to an insulating state as, 
for example, the coupling between the qubits and their 
resonators is increasecP"^", similar to the superfluid- 



Mott insulator quantum phase transition which has 
been observed in ultracold atom systems^ More 
exotic phenomena including analogs of the fractionaP" 
and anomalous^" quantum Hall effects and of Majorana 
physics^ have also been considered. However, while pre- 
liminary steps have been taken to build such a CQED lat- 
tice system^, both establishing and probing its expected 
many-body states remain major experimental challenges. 

The simplest method of probing a microwave circuit is 
to measure transmission between two of its ports. How- 
ever, in the case of a CQED lattice such a measurement 
gives only limited insight into the detailed behavior of 
photons in the interior. Additional information could 
be obtained by measuring transmission while locally per- 
turbing the interior of the sample, as has been done to 
image the coherent flow of electrons in two-dimensional 




Figure 1. Scanning transmon qubit. The qubit (red) is com- 
posed of two 40 x 500 um aluminum islands separated by a 
60 um gap and connected by a thin wire containing a Joseph- 
son tunnel barrier. In the qubit discussed here, the tunnel 
barrier is formed by two aluminum oxide barriers in the two 
arms of a loop in the middle of the wire. The qubit is pictured 
11 um above a portion of a CPWR. The resonator's center pin 
(yellow) is 21 um wide and is separated from the surrounding 
ground plane (blue) by 12 um gaps. The ends of the resonator 
are formed by gaps in the center pin like the one shown in the 
figure. The resonator has a length of 7, 872 um which cor- 
responds to a half-wave resonance of 7.6 GHz in the absence 
of the qubit. The resonator center pin has no curves in over 
half its length, so that the qubit remains aligned as it scans 
over this part of the resonator. To make them more visible, 
the qubit and resonator have been thickened, but all other 
dimensions are to scale with the samples discussed here. 
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electron gas systems^"^]. This local perturbation re- 
quires a new scanning probe tool, such as the one demon- 
strated here. Besides just perturbing the lattice, a scan- 
ning probe qubit can be used to measure the photon num- 
ber of individual lattice sites following a protocol used to 
measure photons in non-scannable cavities^. One bene- 
fit of a scanning qubit in this case is its ability to mea- 
sure the photon number of interior lattice sites. Measure- 
ments of outer resonators, those most easily accessed by 
a measurement circuit fabricated on the same chip as the 
lattice, would be difficult to interpret due to edge effects. 

While our primary focus is on using the scanning 
qubit in the context of CQED lattice-based quantum 
simulation, we note that scanning qubits have been 
studied previously for other applications. Scanning 
nitrogen-vacancy (NV) center qubits in diamond have 
been demonstrated to be sensitive local probes of mag- 
netic nekP^^. A diamond NV center has also been 
coupled to a scanning photonic crystal cavity^. While 
this experiment did not reach the strong coupling regime, 
the scanning cavity enhanced spontaneous emission from 
the NV center and allowed its position to be determined 
with greater spatial resolution. 

The scanning qubit described here (Fig. [I]) is a trans- 
ition design consisting of two aluminum islands connected 
by a thin aluminum wire interrupted by an aluminum 
oxide tunnel barrier^. The tunnel barrier provides a 
large nonlinear inductance, which together with the ca- 
pacitance between the two islands, makes the transmon 
behave as a nonlinear LC oscillator whose lowest two en- 
ergy states can be used as a qubit. The transmon design 
is well suited for scanning because it couples to CPWRs 
capacitively and requires no physical connections. The 
qubit chip was mounted face down to a cryogenic three- 
axis positioning stage and positioned over a separate chip 
containing a niobium CPWR with a half-wave resonance 
at 7.6 GHz. In order to avoid direct contact between the 
resonator and the qubit, pads of photoresist 7/im thick 
were deposited on the corners of the qubit chip. The sam- 
ple holder was mounted to a dilution refrigerator which 
operated at temperatures < 35 mK. 

The main result presented here is the measurement of 
the strength g of the coupling between the resonator and 
the qubit as a function of qubit position. Following ref. 
l2"4l the Hamiltonian H describing the coupled resonator- 
qubit system can be approximated by 

H = hv r ^a f a + + ^~^ z + ^ (aa + + a f <7~) (1) 

with v r and v q the resonator and qubit frequencies re- 
spectively. In this expression, a, are the creation and 
annihilation operators associated with photons in the res- 
onator and <7 + , <t~, and a z are the Pauli spin matrices 
associated with the qubit when treated as a two-level sys- 
tem. On resonance (y q — v r ), the first two excited states 
of the system are (|0 f) ± |1 4-))/v2 with corresponding 



energies hv r ± hg above that of the ground state |0 4-) 
where \nq) is the state with n photons in the resonator 
and the qubit in state q with I (t) representing the qubit 
ground (excited) state. When driven with a microwave 
excitation, transitions to each of these excited states are 
allowed, resulting in two peaks in the low power transmis- 
sion spectrum. These peaks are separated in frequency 
by 2g, called the vacuum Rabi splitting. 

The frequency v r of the resonator depends on its ca- 
pacitance to its ground plane which is decreased when 
the qubit chip is brought into close proximity. In or- 
der to ensure that resonance was possible at every qubit 
position, the qubit was fabricated with a pair of tunnel 
barriers integrated into a loop in place of a single tun- 
nel barrier. By varying the flux through this loop with 
a magnet coil incorporated into the positioner, the qubit 
frequency v q could be varied from a maximum value of 
12.1 GHz to close to zercB^. We note that, although here 
the flux loop's only purpose is to make the qubit energy 
tunable, such a loop can also be operated as a sensitive 
local magnetometer in a scanning SQUID microscope^. 

Fig. [2] shows the transmission spectra of the resonator 
for a sequence of regularly spaced qubit positions along 
the y axis perpendicular to the resonator. At each posi- 
tion, the current through the magnet coil was adjusted 
to bring the qubit into resonance, where the single trans- 
mission peak of the resonator was transformed into two 
peaks of equal height, clearly demonstrating strong cou- 
pling between the scanning qubit and the resonator. The 
position scan shows two regions of large peak separa- 
tion symmetric about a position with nearly no peak 
separation which we set as the origin. In coupling to 
the resonator the transmon behaves as a dipole antenna. 
Because the two islands of the qubit are identical, by 
symmetry no coupling is expected when the qubit is cen- 
tered above the resonator at y — 0. The points of maxi- 
mum peak separation occur at y ~ ±50 um where one of 
the two islands is centered over the resonator. At these 
points, the observed coupling strength g ps 140 MHz was 
well into the strong coupling regime g > k,T^ 1 where 
the qubit relaxation time T% — 2.6 us was determined by 
time-domain measurements (see Supplemental Methods) 
and the photon escape rate k = 13 MHz was determined 
from the resonator linewidth. The photon escape rate 
was chosen to be large in order to increase the rate of 
data acquisition. 

Scans of resonant transmission versus y position like 
Fig. [2] were repeated at five positions along the length of 
the resonator (the x direction) with a spacing of 600 um. 
The coupling strengths g extracted from fits to the trans- 
mission spectrum at each qubit position are plotted in 
Fig. [3] The coupling strength increases as the qubit 
moves from the voltage node at the center of the res- 
onator to the antinode at its end but exhibits the same 
shape for its y dependence at each x position. 

We now consider the data in Fig. [3] more quantita- 
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Figure 2. Resonant transmission spectra. Transmission spectra taken with the qubit in resonance with the resonator are shown 
for a series of qubit positions. The vertical axis represents the difference between the frequency of the applied microwave drive 
and the resonator frequency. The shifted peak locations can also be interpreted as values of the coupling strength g which 
is equal to the peak separation. We take the y origin to be the point of smallest peak separation which we interpret as the 
position where the qubit is centered over the resonator. The suppressed transmission at positions y — ±30 and ±125 um is due 
to coupling between the resonator and a parasitic resonance in the metal frame of the qubit chip (see Supplemental Discussion) . 
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Figure 3. Coupling strength g versus qubit position. Traces of g versus qubit y position are shown for five qubit x positions 
spaced 600 um apart from each other. The traces correspond to scans like the one shown in Fig. [2] Each value of g was 
determined by fitting several transmission spectra taken at values of magnetic flux for which the qubit frequency was close to 
that of the resonator. The sign of Ax was chosen so that with increasing Ax the qubit moved from the electric field node at 
the center of the resonator towards the electric field antinode at its end. 
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Figure 4. Quantitative analysis of measured coupling 
strength, a, The maximum coupling strength g of each trace 
in Fig. [3] is plotted versus x position along with a fit to the 
expected sinusoidal dependence. The x origin represents the 
midpoint of the resonator. The offset of the data points from 
the resonator midpoint was determined by the fit (see Sup- 
plemental Methods). The other fit parameter, the maximum 
coupling strength the end of the resonator, was found to be 
185 MHz. b, The largest trace of g vs. y in Fig. [3] is replot- 
ted along with a fit to the form expected from finite element 
modeling of the qubit-resonator system's capacitance matrix. 
The fitting function uses the resonator frequency u r , the sys- 
tem geometry, and the qubit x position determined in panel 
a as fixed inputs and treats the qubit height z, found to be 
11 um, as its only free parameter. 



tively. The coupling strength g is proportional to the 
electric field at the position of the qubit when a single 
photon is present in the resonator. Panel a of Fig. [4] 
shows the maximum value of g for each x position along 
with a fit to the expected sinusoidal dependence of the 
electric field strength along the x direction. The fit pro- 
vides an absolute reference for the relative x positions 
quoted in Fig. [3] With the x position known, we can 
compare the measurement data with the g(x,y,z) ob- 
tained from finite element simulations (see Supplemen- 
tal Methods) . Panel b of Fig. [4] shows the values of 
g observed for x = 3, 330 um from the resonator mid- 



point along with the simulation results which show good 
agreement. The simulation height z = Hum is some- 
what larger than the 7 pm thickness of the photoresist 
pads on the corners of the qubit chip but corresponds to 
a misalignment between the qubit and resonator chips of 
~ 0.1° over the 2 mm from the edge of the qubit chip to 
the qubit's location at the chip center. 

In conclusion, we have observed strong coupling be- 
tween a scanning transmon qubit and a CPWR. Because 
strong coupling can be reached, anything possible with 
stationary qubits now becomes possible with a moveable 
probe opening the door for a wide array of applications 
in the field of CQED. Such a scanning qubit makes pos- 
sible quantum measurements of superconducting circuits 
with spatial resolution. In addition to the scanning mea- 
surements of a lattice CQED system discussed earlier, 
we note that the system studied here demonstrates in 
situ tuning of the coupling strength g which is often a 
desirable capability experimentally. While the scanning 
qubit's coupling can not be tuned on the timescale of the 
coherence time like some previously demonstrated circuit 
designs^ ^1, it does not require flux or current biasing of 
the system. For example, for the lattice CQED system 
mentioned above, a lattice of resonators on one chip could 
be coupled to a lattice of qubits on a second chip allow- 
ing the coupling between each resonator-qubit pair to be 
tuned together as one chip is scanned over the other. An 
array of qubits could also be used to measure the statis- 
tics of qubit coherence by scanning the qubits one by one 
across a single measurement resonator. 

METHODS 

The qubit was fabricated using electron beam lithog- 
raphy and double-angle shadow evaporation with con- 
trolled oxidation of 30 and 100 nm layers of aluminum 
onto a 4 x 4mm sapphire chip. The 0.5 x 1.0mm crash- 
pads on the corners of the chip were made with pho- 
tolithography of SU-8 2005 photoresist. The resonator 
was defined by photolithography and acid etch (H2O, 
HF, and HN0 3 in a 7.5:4:1 ratio) of a 200 nm film of 
niobium on a 14 x 14 mm sapphire chip. 

The qubit chip was glued with methyl methacrylate to 
the tip of a highly conductive copper rod mounted to the 
cryogenic positioning stage (Attocube ANPx340/RES, 
ANPzlOl/RES). The resonator chip was mounted to a 
copper patterned circuit board with silver paste and alu- 
minum wire bonds which connected the input and output 
transmission lines to coaxial lines. Wire bonds were only 
placed around the edge of the chip outside the footprint 
of the qubit chip. The wiring scheme of the coaxial lines 
was the same as that described in ref. [55] 

All values of qubit position were determined by po- 
tentiometric measurements of resistive position encoders 
integrated into the positioning stage. Individual position 



5 



readings had an uncertainty of 0.4 una. and overall the po- 
sition readings drifted by 1.8 um per 100 um. A typical 
movement of the positioning stage heated the refrigerator 
from its base temperature of < 15 mK to over 85 mK. In 
order to reduce measurement time, most measurements 
were taken with the refrigerator in the range between 25 
and 35 mK which took only a couple minutes to reach 
after moving the stage. 
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SUPPLEMENTARY METHODS 
Fitting functions 

Non-linear least-squares fitting routines were used to 
determine the coupling strength g from transmission data 
and to produce the curves shown in Fig. [4] Here we 
briefly describe the functions used in each case. 



Transmission measurements 

Transmission measurements were made in the low 
power limit for which the rate of photons entering the 
resonator was less than the escape rate, so that the res- 
onator occupancy was less than one photon on average. 
In this case, the transmission spectrum contains peaks 
at frequencies v± corresponding to transitions from the 
ground state |0 J.) to the states |+) and |— ) which pos- 
sess a |1 I) component and are located at energies hu + 
and hv_ above the ground state. The states and their 
frequencies are found by diagonalizing the Hamiltonian 
given in equation 0: 




and 



v± = 



2 



A 2 



(SI) 



with A = v q — v r . The peak amplitudes are pro- 
portional to the probabilities w± = |(1 J, l±)| 2 of a 
photon being measured in the states |±). The peak 
linewidths j± are equal to the decay rates of the qubit 
(Tj -1 ) and the photon (k) weighted by the probability of 
measuring a qubit excitation and a photon respectively: 
7± = w±n + (1 — w±)T^ . We take transmission peaks 
as following a lorentzian lineshape and use the following 
form to fit the resonator transmission: 



S 21 {v)=B + A 



(S2) 



where A is the overall amplitude accounting for all at- 
tenuation and amplification in the measurement circuit, 
B is the background of the detector, and l(v, 7) is the 
complex lorentzian centered at i>q with width 7: 



7/2 



When fitting for g, the parameters A, B, k, v qi u r , and g 
were allowed to vary, while T\ was held fixed to the value 
obtained from coherence measurements. Fig. |S1| shows 
the result of fitting one of the transmission spectra from 
Fig. [2j Also shown is a plot of transmission versus flux 
as the qubit passes through resonance. For the coupling 
strength values shown in Fig. [3j similar flux scans were 
taken and the plotted values of coupling strength were 



obtained by averaging the coupling strength values ob- 
tained from fitting the transmission at each flux value. 



Coupling strength versus x 

The voltage profiles of the modes of a CPWR with 
open boundary conditions are sinusoidal along the length 
of the resonator with antinodes at its ends. The coupling 
strength is proportional to the resonator voltage with one 
photon present and so should follow this sinusoid. The 
maximum coupling strength at each x position shown in 
panel a of Fig. [4] was fit to the sinusoidal form: 



g (Ax) = g n 



sin 7T 



Aa 



x 



(S3) 



where l r = 7, 872 pm is the resonator length and <7 ma x 
and Xq were the fitting parameters. Here Ax is the set 
of displacements in x from the first x position (i.e. the 
values are Oyrm; 600um; l,200um; etc.). In panel a of 
Fig. [4] the measured coupling strengths and the fit are 
plotted versus x — Ax + xq. 



Coupling strength versus y 

In order to perform the fit of coupling strength versus 
y shown in panel b of Fig. |4j the following expression for 
the coupling strength given in ref. [53|was used: 



g(x,y,z) = 2 



2Z r 



m (x) v r (3 (y, z) n al (y, z, v r 
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Figure SI. Fitting resonator transmission, a, The transmission at y — 146 pm in Fig. [2]is plotted in arbitrary units along with 
a fit to equation (S2|. The fit coefficients are A = 57, B = 0.2, g = 20 MHz, k = 14 MHz, v r = 8.342 GHz, and v q = 8.339 GHz. 
b, Transmission spectra recorded at the same position as panel a are plotted versus flux in arbitrary units as the qubit passes 
through resonance. Averaging the coefficients from fits to the transmission at each flux value gives A = 59 ±3, B = 0.30 ±0.05, 
g = 20.4 ± 0.3 MHz, k = 12.8 ± 0.8 MHz, and v r = 8.342 ± 0.001 GHz where the errors represent the standard deviation of the 
coefficients. 



with the characteristic line impedance Z c taken to be 
50 fl and m(x) the sinusoidal mode shape factor given in 
equation ( |S3[ ). In order to describe the voltage division 
factor /? and the transmon matrix element noi, we first 
define Cjk to be the capacitance between components 
j and k and label the components of the system with 
a and b for the two islands of the transmon, p for the 
resonator center pin, and g for all other pieces of metal 
(the two ground planes and the metal frame on the qubit 
chip). The voltage division factor j3 gives the fraction of 
the voltage drop from the resonator center pin to ground 
that falls across the two islands of the qubit. It can be 
written in terms of capacitance coefficients as 

P _ \C a pCbg — Ct,pC a g\ 

Cab (C a E + Cf,s) + C a sCbT, 

where C x s = C xp + C xg . We find the matrix element uqi 
by numerically diagonalizing the transmon Hamiltonian 
given in ref. [2H 

H = AE c h 2 - Ej cos (p 

to finds it eigenstates and eigenenergies and then evalu- 
ating n i = (0|n|l), where |0) and |1) are the eigenstates 
with the two lowest energies, E and E\. The charging 
energy Ec — e 2 /2Cs used in the calculation was calcu- 
lated using the total capacitance given by 

The qubit frequency v q is given by (Ei—E a )/h and is thus 
a function of Ec and Ej. In calculating noi, v q (Ec, Ej) 
was numerically inverted to solve for Ej(Ec,v q ) with 



v q set equal to v r since measurements of the coupling 
strength were made with the qubit close to the res- 
onator's frequency. 

In order to produce the fit shown in panel b of Fig. 
[4] the coupling strength g(x,y,z) was calculated using 
the known values of Z c and v ri the value of x obtained 
from the fit in panel a of Fig. |4j and the values of the 
capacitances Cjk found by finite element analysis for a 
grid of y and z values with 1 pm spacing. The measured 
coupling strength versus y was fit to the g(y,z) found 
by interpolating between the y and z grid points with z 
as the only free parameter. The finite element simula- 
tion was then repeated with the fitted value of z in order 
to produce the curve shown in panel b of Fig. [4] We 
note that at the fitted value of z = 11.0 pm the charg- 
ing energy Ec — 388 MHz is similar to values used in 
other CQED experiments and corresponds to a ratio of 
Ej/Ec = 59, within the transmon regime where the off- 
set charge across the transmon islands (not included in 
the Hamiltonian given above) may be ignored. 

By the use of alignment marks on the resonator and 
qubit chips, it was possible to confirm that the misalign- 
ment between the two chips was 3° ± 1°. A misalignment 
of 3° was used in the finite element calculations for the 
capacitance coefficients. Using a misalignment of 2° (4°) 
instead gave a the fitted height z of 11.1pm (10.7pm). 

Coherence measurements 

Qubit coherence times (T\ = 2.6 ± 0.3 pis, T 2 * = 1.0 ± 
0.2 ps) were obtained using the techniques described in 
ref. 1291 The measurements were made during the same 
cooldown and at the same x position as the data shown 
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in Fig. [2] For technical reasons, the measurements were 
made immediately after the refrigerator was warmed up 
to 20 K and then cooled back down to its base temper- 
ature. The coherence measurements were performed at 
y = -113 mn (g = 31 MHz) in Fig. [2] with the qubit 
frequency detuned 700 MHz below the resonator. 



SUPPLEMENTARY DISCUSSION 

Here we provide some additional measurements and 
analysis related to scanning the qubit over the resonator. 
In Fig. [2] the transmission is suppressed near y = ±30 
and ±125 um due to the resonator's coupling to spurious 
modes. Because the behavior of these modes is symmet- 
ric in qubit position, we believe them to be caused by res- 
onances between the resonator chip and a layer of metal 
on the qubit chip that was patterned symmetrically. The 
extra metal on the qubit chip was deposited for techni- 
cal reasons and is not needed for the functioning of the 
qubit. By redesigning the qubit chip or resonator chip, 
these modes could be eliminated. 

In panel a of Fig. |S2[ the data from Fig. [2] is replot- 
ted after removing the background and normalizing the 
transmission peaks to unity in order to make the peaks 
near y — ±30 and ±125 um more visible. In Fig. [2] 
the origin of the frequency axis was set to the location 
of the high power transmission peak and the qubit fre- 
quency was tuned to produce two peaks of equal height 
in transmission. For some positions near where the res- 
onator coupled to the spurious modes, these conditions 
produced two peaks not centered around zero. In panel 
a of Fig. [S2] the frequency axis at each position has 
been shifted to center the peaks around zero, so that the 
transmission peaks at neighboring positions can be more 
easily compared. 

The spurious modes appeared as lower and broader 
peaks in transmission at frequencies that varied with po- 
sition and did not vary with magnetic flux. These modes 
were always present but only affected the measurement of 
the resonator-qubit system when their frequencies were 
close to the resonator frequency. When the frequency 
of one of these modes was close to the resonator fre- 
quency, the coupling between the spurious mode and the 
resonator resulted in two modes with excitations partially 
of the resonator and partially of the spurious mode. The 
narrower peak of the resulting two peaks in transmission 
was chosen to be the resonator peak for the purpose of 
coupling to the qubit. Panel b of Fig. S2 shows the fre- 



quency of the chosen peak for each position in Fig. [2] 
Jumps in the resonator frequency due to avoided cross- 
ings with spurious modes are visible at the y positions 
with low transmission in Fig. [2] 

The scan shown in Fig. [2] was taken on a separate 
cooldown from the scans shown in Fig. [3] The same 
resonator and qubit samples were used for both sets of 



measurements, but the sample stage was disassembled in 
between the cooldowns. During the cooldown in which 
data in Fig. [2] was taken, the qubit's x position was 
not varied, so the absolute x position of the data is not 
known. However, using the maximum value of g from 
Fig. [2] and the curve shown in panel a of Fig. [4] to 
calibrate the x position, one finds the data in Fig. [2] was 
taken at x = 2, 116 pm. 

Changing the position of the qubit also affected the 
threading of flux through the qubit's SQUID loop. Once 
the SQUID loop was moved away from the gaps in the 
coplanar waveguide and positioned above the supercon- 
ducting ground plane, the amount of flux produced by 
the magnet coil required to tune qubit into resonance 
increased rapidly because the Meissner effect screened 
the magnetic field away from the superconducting ground 
plane. In panel c of Fig. |S2] the coil magnetic flux that 
brought the qubit into resonance with the resonator is 
plotted for each position of Fig. [2] The impact of the 
magnetic field screening could be greatly reduced by fab- 
ricating holes in the resonator's ground plane. 

The strong dependence of the resonant magnetic flux 
on the qubit position as well as the steepness of the slope 
of qubit frequency versus magnetic flux (maximum qubit 
frequency u q ^ max ~ 12.1 GHz) necessitated careful flux 
scanning at each qubit position in order to locate reso- 
nance. Searching for resonance by monitoring the trans- 
mission spectrum for an avoided crossing feature like the 
one shown in panel b of Fig. SI would have required a 



long measurement time at each qubit position. Instead, 
only the low power transmission at the frequency v r of 
the high power transmission peak was monitored as the 
flux was swept. 

When \v r — v q \ ^> g, one of the two mode frequencies 
v± given in equation (SI I differs from v r by ~ g 2 / {v r — v q ) 



and is associated with a large peak in transmission. For 
most qubit frequencies this frequency shift is small com- 
pared to the resonator linewidth n = 13 MHz and so 



transmission at v r is high. However, when v q 



the 



mode frequencies v± are shifted from v r by g > k and 



transmission at v r is low. Panel a of Fig. S3 illustrates 
this behavior by showing showing transmission at v r ver- 
sus flux and input power. Regular dips in transmission 
occur at low power where the qubit passes through res- 
onance. Panel b of Fig. S3 plots just the low power 



transmission versus flux and shows that resonance can 
be easily identified by monitoring transmission at just 
one frequency value. In practice, a scan like that shown 
in panel b was taken at each position to identify the reso- 
nant flux range, and then a scan like that shown in panel 



b of SI was taken over this flux range to obtain trans- 
mission spectra to fit for g. 

Additional features are present in the crossover from 
the low power region to the high power region of the 
transmission. Panel c of IS II shows a finer scan of trans- 
mission versus power and flux at a qubit position close 
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Figure S2. Additional information about Fig. [2] a, The transmission spectra with the qubit tuned into resonance with the 
resonator is plotted for a series of qubit positions. At each position, the background is removed and the transmission is scaled 
so that the maximum transmission is unity. The frequency of each spectrum is offset so that the two peaks are centered around 
zero, b, The resonator frequency is plotted for each qubit position of Fig. [2] The frequency axis of each spectrum in Fig. [2] 
was offset by the frequency plotted here, c, The total flux (in arbitrary units) sourced by the magnet coil to bring the qubit 
into resonance is plotted versus qubit position for the scan shown in Fig. [2] 



to that of the scan in panel a. These features are likely 
related to higher level qubit transitions coming into res- 
onance with the resonator, though additional analysis is 
needed. 

In the main text, results have been presented for the 
coupling of the qubit to the resonator as a function of lat- 
eral position (x and y). Measurements of the coupling's 
dependence on the qubit 's vertical displacement z from 
the resonator were not possible due to misalignment be- 
tween the resonator and qubit chips. Evidence of this 
misalignment is visible in panel a of Fig. |S3| which plots 



the resonator frequency versus the positioner's z read- 
ing denoted by z p . The origin of z p was chosen to be the 
point at which the positioner could no longer advance due 
to contact with the resonator chip. Above z p = 40 um, 
the resonator frequency is shifted to higher values as the 
qubit chip is brought closer as expected due to the modi- 
fication of the resonator's effective dielectric constant by 
the qubit's presence. Below 40 pm, the resonator fre- 
quency's dependence on z p weakens and disappears even 
as the positioner continues to move. We interpret this 
behavior as the qubit chip coming into first partial con- 
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Figure S3. Transmission at the resonator frequency versus input power and magnetic flux, a, The transmission at v r is plotted 
versus input microwave power on a log scale and magnetic flux threading the qubit loop. The flux axis has been scaled by 
the observed flux period 0q. The power axis is uncalibrated, but the cross-over near -15 on resonance occurs as the photon 
occupation of the resonator increases from 0.1 to 10. The scan was taken at the same x position as Fig. [2]with the qubit close 
to the center pin and g = 29 MHz. b, The average of the transmission for powers less than -19 in panel a is plotted versus flux, 
c, Another scan of transmission versus power and flux like that in panel a is plotted for a flux region close to resonance. The 
flux axis is uncalibrated. The scan was taken at a position with g = 21 MHz, close to the position used for the scan shown in 
panel a . 



tact with the resonator and then nearly full contact as 
compliance in the sample holder allow the two chips to 
align. We attribute the relatively small magnitude of the 
discrepancy of the qubit height obtained by the fit shown 
in Fig. [4] from the height of the photoresist pads to this 
compliance in the sample holder. 

We note that the shift of the resonator frequency due to 
the qubit chip shown in panel a of Fig. |S4| demonstrates 
the possibility of using a scannable chip to produce a 
defect in a CQED lattice by shifting the frequency of 
one of the resonators within the lattice. For scanning 
experiments where no defect is desired, the shift shown 
in panel a could be made unimportant by scanning with 
a chip larger than the entire CQED lattice. In this case, 
the resonator at each lattice site would receive the same 
shift. 

The resonator frequency's dependence on the qubit 



chip height shown in panel a of Fig. |S4] allowed the res- 
onator to be used to measure the qubit height. Panel b 



of Fig. S4 plots noise spectra of the transmitted phase at 
the resonator frequency when the qubit chip is at z p = 
and 43 um. We interpret the phase fluctuations present 
when the qubit chip is hovering above the resonator and 
not present when the qubit chip is in hard contact as 
being due to motion of the qubit chip relative to the 
resonator chip. Using the resonator's phase versus fre- 
quency curve to convert the phase noise into an effective 
resonator frequency noise and then the curve shown in 
panel a of Fig. |S4]to convert frequency into position, we 
obtain the position noise spectrum shown in panel c of 
Fig. 



S4 



Repeating this procedure at another z p posi- 
tion and using the amplitude noise instead of the phase 
noise resulted in vibration spectra with similar features 
at similar magnitudes, confirming the interpretation of 



11 




CD 
"O 
=3 



-Q 
> 



20 40 60 80 100 120 
z positioner reading (urn) 

"T 




0.01 



100 200 300 

Frequency (Hz) 




600 800 
Frequency (Hz) 



1000 



1200 



1400 



Figure S4. The effects of retracting the qubit from the resonator chip, a, The resonator frequency is plotted versus the reading 
z p of the z positioner. At z v — 0, the qubit chip is in hard contact with the resonator chip and can not be advanced further, 
b, The noise spectrum of the transmitted phase at the resonator frequency v r is plotted for the positions z v = and 43 p. c, 
The noise spectrum of the qubit-resonator displacement is plotted for z p = 43 pm using the data from panel b. The slope of 
the curve shown in panel a was used to convert the phase data into displacement. 



the features as being due to vibration of the qubit chip. 
The spectrum shown in panel c of Fig. |S4| is typical for 
the mechanical response of cryogenic positioners such as 
those used in this experiment. The motion of the re- 
frigerator base plate inferred from the spectrum shown 
in panel c of Fig. |S4| agreed with measurements made 
with an accelerometer of another refrigerator that was 
the same model as that used for the measurements pre- 
sented here. Because all measurements of the qubit were 
made with the resonator and qubit chips in hard con- 
tact, no vibration isolation elements were included in the 
sample holder. In order to measure the dependence of 



the qubit-resonator coupling on height, such vibration 
isolation would need to be considered in addition to the 
alignment of the qubit and resonator chips. 
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